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Abstract 

We present a theory of resonances for a class of non-autonomous Hamiltonians to treat the 
structural instability of spatially localized and time-periodic solutions associated with an 
unperturbed autonomous Hamiltonian. The mechanism of instability is radiative decay, 
due to resonant coupling of the discrete modes to the continuum modes by the time- 
dependent perturbation. This results in a slow transfer of energy from the discrete modes 
to the continuum. The rate of decay of solutions is slow and hence the decaying bound 
states can be viewed as metastable. The ideas are closely related to the authors' work on 
(i) a time dependent approach to the instability of eigenvalues embedded in the continuous 
spectra, and (ii) resonances, radiation damping and instability in Hamiltonian nonlinear 
wave equations. The theory is applied to a general class of Schrodinger equations. The 
phenomenon of ionization may be viewed as a resonance problem of the type we consider 
and we apply our theory to find the rate of ionization, spectral line shift and local decay 
estimates for such Hamiltonians. 
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1. Introduction 



There are many dynamical systems of interest for which the dynamics can be viewed as given 
by a Schrodinger type equation generated by an operator H(t) = H Q + W(t). Here, H is a 
time-independent (autonomous) part, assumed to be self-adjoint on a complex Hilbert space 
Ti, and W(t) denotes time-dependent (non-autonomous) perturbation. We consider the initial 
value problem for the Schrodinger equation for <fi : t \— > <p(t) G Ti: 



with initial condition 0(0) = 0o m a suitable subspace of Ti. An important special case, studied 
in section 3, is: 



where H = —A + V(x), with V(x) real-valued and W(t, x) is a smooth function which is 
spatially localized and time-periodic or quasi-periodic. Equations of this type play an important 
role in tunneling and ionization physics, and are relevant to an understanding of many nonlinear 
problems. This is discussed later in this section. 

In this paper, we consider the situation where the unperturbed equation, corresponding to 
the case W = 0, has a time-periodic solution in Tl, e.g. a bound state. We study the structural 
instability of the bound state due to a perturbation W(t), which is assumed to be periodic. 
The analysis is carried out for the special case W(t) = cos(/xt) (3, with (3* = (3; see section 2. 
The case of more general perturbations W(t), which are periodic or quasi-periodic, as well as 
the case when the unperturbed operator Hq has more than bound state can be treated by our 
approach. We indicate briefly in section 4 the modifications required in the analysis. 

Viewed in terms of the coordinates of the unperturbed dynamical system, the perturbation 
W(t) has the effect of coupling the bound state to the continuum modes. The mechanism of 
instability is a slow transfer of energy from the bound state part of the solution to the continuum 
modes and an associated decay due to radiation of energy to infinity. We develop a theory which 
yields a detailed picture of the intermediate and long time behavior of solutions to the initial 
value problem. Our main results are presented in Theorems 2.1 and 2.2. A consequence is that 
the solution of the initial value problem for the perturbed dynamical system with initial data 
given by the time-periodic bound state of the unperturbed problem, is characterized by transient 
exponential decay ~ e _r '*' where < T = 0(\ \ \W\ | | 2 ), for |||W||| small, and by algebraic decay 
as \t\ — > oo. We derive explicit formulae for the lifetime, t = T^ 1 (eqn. Q2.9|) ) and the spectral 
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(H + W(t))4>, 






line shift (eqn. (|4.21|) ). 
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The approach we take is closely related to our work on quantum resonances |57| , [38) and on 
the radiation damping of localized states for nonlinear wave equations |39[]. We now illustrate 
the main idea in the current context. Consider a linear Schrodinger equation: 



id t <j>(t,x) = (-A + V(x))(f)(t,x) + cos(pt) /3(x)4>(t,x). (1.3) 

We assume that —A + V(x) has exactly one eigenvalue at A < 0, with corresponding eigen- 
function ipo(x). Therefore, the unperturbed (fl = 0) problem has a spatially localized and 
time-periodic solution, ilj (x)e~ tXot . We assume that the continuous spectrum of H is equal 
to the real nonnegative half-line. The perturbation W(x,t) = (3(x) cos(/x£) is assumed to be 
real- valued, to decay sufficiently rapidly as \x\ — ► oo, and to be small in an appropriate norm; 
see section 3. 

We now ask the following question: Does the bound state solution of the unperturbed problem 
persist under small perturbations of this type? 

As indicated above, we find that generic perturbations lead to a slow radiative decay of solu- 
tions. To illustrate the mechanism of decay, we first observe that it is natural in the unperturbed 
case to decompose the solution as: 

<f>(t,x) = a(t)Vo + Mt,x), (ipoAdit)) = 0. (1.4) 

Then, we have that 

id t a(t) = A a(t), 
i<j> d (t,x) = (-A + V(x)) (f> d (t,x), (1.5) 

with initial conditions 

a(0) = (Vo,0(O)), <M0) = Pc0(O) = 0(0) - (Vo,0(O))Vo. (1.6) 

P c is the projection onto the continuous spectral part of the operator H . The equations (|1.5|) 
decouple and we have that all solutions of ( |1.4|) are of the form (|1.4|) with: (i) a(t) = e~* A °* a(0), 
and under suitable hypotheses on V(x), (ii) 4>d(t,x) decays (dispersively) as |t| — > oo, in the 
sense that ||^d(t, -)\\ p — > 0, (p > 2) or for some a > 0, || (x)~ a (f) d (t, ■) \\ L 2 — > 0. 

How is the dynamics effected by a small spatially localized and time-periodic perturbation 
(P 7^ 0)? Substitution of ( |1.4j ) into (|1.3|), leads to a weakly coupled dynamical system for a(t) 
and <pd{t)- This system is derived and studied in sections 4-6. A simplification of the resulting 



3 



system which we shall use to illustrate the main idea is the following: 

id t a(t) = A a(t) + ( f3ip , 4>d{t) ) cos(/it) , 
id t <f>d(t) = -N>d(t) + a(t) cos(fit) (3^o- (1.7) 

The system (|1.7| ) can viewed as an infinite dimensional dynamical system consisting of two 
components: (i) a finite dimensional part, described an "oscillatory", with complex amplitude, 
a(t), coupled to (ii) an infinite dimensional dispersive wave equation on IR n , governing (f>d(t,x). 
It is simple to check that the system ( |1.7|) has the conserved energy: 

£[a,<j> d ](t) = \a(t)\ 2 + J \Mt,x)\ 2 dx = 5[a(O),0 d (O)]. (1.8) 

The first term in ( |1.8|) corresponds to the oscillator energy and the second to the wave energy. 

Intuitively, a spatially localized finite energy disturbance of the rest state should result 
in motions induced in the oscillator as well as waves propagating "outward" and eventually 
escaping any compact set (a characteristic of dispersive waves). Due to the coupling (/3 ^ 0) of 
the oscillator and waves, we expect the energy of the oscillator, |a(t)| 2 , to decay toward zero as 
time progresses. This energy is transferred to the continuum modes, while the total energy of 
the system, £[a, u], remains constant in time. This is what is proved under suitable hypotheses 
in Theorems 2.1 and 2.2. 

The way in which we establish energy transfer from the discrete to the continuum modes 
is by showing that one can transform a system of the type (|1.7| ) to a normal form, in which 
internal damping, which drives the energy transfer, is made explicit. This damping is due to a 
resonant coupling of the discrete oscillator mode with the continuum modes. We now explain 
the origin of this effect. 

Viewing (3 as small, we set a(t) = e~ lXot A(t), where A(t) is "slowly varying", (dfA{t) = 
O(P)). In terms of the amplitude, A(t), the system (|1.7f ) becomes: 

d t A(t) = -i ( , M*) ) cos ^ eiAo *> (I- 9 ) 
id t <j) d = -A cf) d + A(t) e~ iXot cos(^)/3 % (1.10) 

Since the variation of A(t) is slow, the source term in ( |1.10| ) is essentially quasi-periodic in time. 
One of its frequencies is Xo + fi, which we assume to be strictly positive, i.e. A + /U G o cont {— A). 
That is, the non-autonomous parametric forcing of frequency fi induces a resonant forcing of 
the wave- field 4>d- The character of solutions to ( |i~9"l - |1.10|) can b e deduced by first solving ( |1.10| ) 
for (f)d and then substituting the result into ( |1.9| ). The result is a reduction to an ordinary 
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differential equation for A(t). This is reminiscent of the center manifold approach, commonly 
used in dissipative systems (see, for example, ||, [[ll]]) and recently applied in a Hamiltonian 



context |16 |. 

The asymptotic behavior of the ordinary differential equation governing A(t) can be deduced 
by a asymptotic evaluation of the right hand side of fll.9|), considering carefully the contribution 
to <pd coming from a neighborhood of the resonant energy, Aq + /i. 



-i ( PiJo , <f> d (t) ) cos /it e iXot 



e iA(t-s) e -i(\ 0+IJ ,)s A ^ ds 



3 i(Ao+//)t 



+ nonresonant terms 



= — — ( /3i[)o , (— A — Ao — n — iO) 1 (3ipo^j A(t) + nonresonant terms 
Substitution of this result into the equation for A{t) yields: 

d t A(t) = (-r + iA + p(t) ) A(t) t > 

where 



;i.i2) 
;i.i3) 



r = vr ( P4> , 5(-A - Ao - /i) (3% ) . 

A is a real constant (the Lamb shift) and 0{f3 2 ) and p(t) is a bounded oscillatory function of t 
which is 0{(3 2 ). 

The right hand side of (|1.13|) is defined using the spectral theorem for —A on L 2 (IR n ). In 
terms of the Fourier transform, JF[^](^) = J e~ l ^' x z(x) dx, we find: 



Wo 



A + A* 



+ 



T\^o 



A + /i 



(1.14) 



The constant T is always non-negative. If it is strictly positive, then (|1.12 ) is a damped oscillator, 
with damping coefficient T resulting from coupling of the discrete and continuum modes. In our 
general theorems, Theorems 2.1 and 2.2, we have hypothesis (H6), that the analogous general 
expression for T is strictly positive. This holds generically in the space of perturbations W 
under consideration. 

Emergence of irreversible behavior in a time-reversible system 

Thus, the original dynamical system, ( |1.3|) , which is Hamiltonian and, in particular, time- 
reversible ( invariant under the transformation 0(t, x) h- > (f>(—t, x) ), is equivalent to a dynamical 
system with some apparent time-asymmetry. [] This is related to the e - prescription discussed 
in section 4; see, in particular, Proposition 4.2. At the heart of this apparent paradox is that 
1 This is related to a "paradox" raised in the late nineteenth century by Zermelo and others concerning the 
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the energy of the oscillator, which is coupled to the continuum, is propagated out to infinity 
because the singular limit: 

lime-UM (-A - (A + n) =f ^ (1-15) 

eJ.0 



satisfies a local energy decay estimate as t — > ±00; see (H4). Thus, for t < 0, in ( |1.12|) we 
replace T by — T. 

Normal Forms for coupled oscillator - wave systems 

The damped oscillator above is a normal form equation, which illustrates a more general 
phenomenon. Dynamical systems like can be viewed as a Hamiltonian systems consisting 
of two subsystems: a single (or, in general, multiple) degree of freedom oscillator and an infinite 
dimensional radiating wave field. 

(I) In fll .3D the discrete oscillator frequency does not lie in the continuous spectrum of wave 
frequencies. However, the time periodicity of the potential generates frequencies in the solution 
which resonate with continuum modes. 

(II) In Ift^l , [ 3S|] , quantum resonance type problems are considered in which the discrete oscillator 



frequency coincides with a frequency of continuum modes. Coupling of the two occurs due to a 
time-independent localized potential. 

(Ill) In our study of radiation damping of "breather-like" structures we find, as in (|1.3| ), that 
the discrete oscillator frequency does not lie in the continuous spectrum of wave frequencies. 
Here, the coupling is due to a nonlinear potential which generates frequencies which lie in the 
continuous spectrum. 

In an isolated single degree of freedom Hamiltonian system, the general normal form is: 

A' = z(do + c 21 \A\ 2 + c 32 \A\ 4 + ... + c n+1 , n \A\ 2n + ...) A, (1.16) 
where the constants c n+ i jn are real numbers. 

How is this normal form altered as a result of coupling of the oscillator to an infinite dimensional 
radiative wave equation ? 



In the |37]], [38|, [39| and in the present work, we find the more general dispersive or radiative 



Hamiltonian normal form: 

A' = (k 10 + k 21 \A\ 2 + k 32 \A\ 4 + ... + k n+1 , n \A\ 2n + ...) A, (1.17) 

irreversibility of macroscopically observed behavior in physical systems and the reversibility of the microscopic 
equations of motion. This was addressed in the context of classical systems via a statistical approach in the work 
of Maxwell and Thomson and culminated in the work of Boltzmann. In quantum systems, for example in the 
quantum measurement theory, extensions of these ideas apply. The source of the emergence of irreversibility, is 
the coupling of systems with few degrees of freedom to systems with many. For a discussion, see |EQ| . 
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with complex coefficients; 

kn+l,n — d>n+l,n + 2C n+ i j?1 , (1.18) 

where d n +x >n and c ra+ i jn are real numbers. Non-zero real parts of the normal form coefficients 
arise due to coupling of discrete and continuum modes. In the linear problems: (I) and (II), we 
always have dxo = — r < 0. In (III) we have k±o = ic\$ and, depending on the details of the 
nonlinear interaction, we may have: k^x = ic2i and k^ = — r + ic^ j39fl , &2i = — T + ici\ 



or other possibilities. In each case there's a different expression for T, the first nontrivial real 
contribution to a normal form coefficient, but the essence of its calculation is as described above. 
Generically one has — Y < 0, and it is natural to conjecture that the first non-zero T = <i n „ + i jri „ 
is strictly negative. 

Connection with persistence theory of periodic solutions: 

The condition T > (see also (H6)) can be seen to arise in another manner. Motivated by 
the perturbative approach to constructing periodic solutions of ordinary differential equations, 
it is natural to seek a solution of (11.31) in the form: 



(j)(t) = e" iAo Vo + V>i + ^2 + + ^ + -, (1-19) 
where if)j is formally of order /3 J Substitution of the ansatz ( 1.19 ) into (|1.3|) yields that: 



Vi = Vi + (1-20) 

where 



id t ri = Horii + I e-^ Xo+ ^ /3Vo- (1-21) 
These equations can be solved by setting: 

= e -t(Ao-M)t ^ = e -i(\ +n)t ^6 

For ipi we obtain the equation: 

( Ho - (A - n) ) *° = 

(#0 - (A + /i) ) # 6 = ^0 (1-23) 
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The first equation of ( |1.23| ) has a localized solution because Ao— /i < and lies in the resolvent set 
of Hq. Note however, that Ao + fi > 0, and therefore the second equation can have a localized 
solution only if the projection of f3ipo onto the generalized eigenfunction associated with the 
continuum energy Ao + /i vanishes. This is equivalent to Y = 0, where T is defined by ( |1.13| ); 
See also ( |1.14j ) or, more generally, (H6). Generically, T > 0. The results of this paper show, 
in particular, that this obstruction to solvability in the class of quasiperiodic functions, implies 
radiative decay of solutions of the initial value problem. 

Finally, we discuss various contexts, where problems involving time-dependent potentials 
have been studied. Equations with time dependent potentials arise in mathematical physics 
and are important in applications, e.g. electron microscopy and solid state devices. In quantum 
physics their analysis dates back to the early days of the theory, where the process of ionization 
by external field (light waves) was studied and later nuclear spin resonance and charge transfer 
problems; see e.g. [[19], fTj. The problem of ionization and more generally of excitation of a 
molecule by radiation is also important in laser optics and chemistry. Recent models of chaos 



14 and stochastic resonances also involve time dependent potential models, e.g. ||42||. The 



expression ( |1.13|) (see also (H6)) is the analogue of the Fermi Golden Rule, which arises in the 
study of spontaneous emission, autoionization ]7), [19|, [29|. In the context of the problem of 
ionization of an atom by a time-dependent electric field a heuristic derivation is presented in 
section 42 of [19[ . 

Equations of type ( |1 . 1[ ) and ( |1 . 2j ) are also relevant in the study of coherent solutions of 
nonlinear systems. For example, a periodic or quasiperiodic solution can be viewed as a bound 
state of an equation of this type with a self-consistent (nonlinear) time-dependent potential. This 
point of view was first exploited by Sigal |3"3|] , |2"4| who used the methods of dilation analyticity 
to relate the generic structural instability of quasi-periodic solutions of certain nonlinear systems 
to the generic instability of an embedded eigenvalue in the continuous spectrum of a Floquet 
operator. This approach yields the spectral line shift and the Fermi golden rule, which give the 
location of the resonance to which the embedded eigenvalue has been perturbed. The approach 
we have taken to the non-persistence question is to show by a study of the large time behavior 
of the initial value problem that solutions decay slowly as t — > oo. 

The analysis of time dependent resonant problems in the physics literature is based on formal 
time dependent perturbation theory; see |7|], where Dyson series to first order is used, and 
|30H where Dirac's perturbation theory is applied. This analysis is limited in many respects. In 
19] , pi it is assumed, among other things, that the relevant time scale of the approximation is 



smaller than e 1 where e is the perturbation size, while we know that the lifetime, is of order 
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e~ 2 ; It is also not clear how to get the spectral line shift from this approach. 

Problems involving time dependent potentials have received considerable attention in the 
mathematical physics literature during the last twenty years, and a number of different ap- 
proaches, based on scattering theory techniques, have been developed. In J5| a kind of Floquet 
theory is used to analyze the case of time-dependent potentials decaying like L 1 functions in 
time. The approach has also been used to treat some problems involving time-periodic potentials 
and long range potentials [T7|], JEJ, [T0| . 

Another type of time-dependent problem is the charge-transfer problems. These are problems 
in which the potential W is a sum of terms of the type: W(x — vt). In a charge transfer 
problem is mapped in a similar way onto a three-body type problem which is treated by resolvent 
methods; In || a geometric approach to scattering was used to study the general behavior of 
scattering states with time- dependent potentials. 

For potentials decaying slowly in time which are perturbations of many body hamiltonians, 
which appear in the analysis of N-body long range scattering, phase-space methods were used 
in [35], based on the existence of an asymptotic energy operator in such cases. The AC Stark 



effect was treated in [44] by methods exploiting the special form of the potential. Also related 
to our study is research concerning the wave equation with a time periodic potential or the free 



wave equation in the exterior of a periodically oscillating obstacle; see, for example, H, [4C 
and those cited therein. 

We consider a Hamiltonian, Hq, with one bound state, perturbed by a small and localized 
potential, time-periodic potential of the simple form given in Hypothesis (H5). This assumption 
is made to simplify computations but the general method applies to problems with more than one 
bound state and to more general time - periodic or quasi-periodic perturbations. An application 
currently under investigation involves the structural instability and meta-stability of "breather- 
like" modes of certain linear wave-guide problems related to " M-soliton" solutions of integral 
nonlinear flows |H . 

In this paper we treat a model of the ionization problem [19|], [[?[]. By treating the process of 
ionization as a resonance caused by coupling of the point spectrum to the continuous spectrum, 
we apply the time dependent resonance theory, recently developed in p7H , ]|3~8| to find the large 
time behavior of the solution. In particular, we show that under our conditions, the bound state 
is always disintegrating and find the lifetime, r ~ T -1 , and the spectral line shift, ( 4.21|) ; see 
the remark following Theorem 2.2. We also prove local decay estimates for the solution, which 
are of general interest to nonlinear and other problems. 
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Notations and terminology: 

Throughout this paper we will denote a generic constant by C, D, etc. 
(x) = (1 + \x\ 2 )K 

C(TC) = the space of bounded linear operators on Ti 
Functions of self-adjoint operators are defined via the spectral theorem p8 
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tion. 



2. General formulation and main results. 

Consider the general system 

id t (j>{t) = (H + W(t))^t), (2.1) 

Here, <ft(t) denotes a function of time, t, with values in a complex Hilbert space 7i. 
Hypotheses on H : 

(HI) H is self-adjoint on H and both H and W(t), t G IR 1 , are densely defined on a subspace 
V of U. 

The norm on TL is denoted by || • ||, and the inner product of /, g G H, by (/, g). 
(H2) The spectrum of H is assumed to consist of an absolutely continuous part, u C ont(H ), 
with associated spectral projection P c and a single isolated eigenvalue A with corresponding 
normalized eigenstate, ipo, i-e- 

H ijj = Ao^o, HVoll = 1. (2-2) 

The manner in which we shall measure the decay of solutions is typically in a local decay 
sense, e.g. for the scalar Schrodinger equation we measure local decay using the norms: / i— > 
||(x) _fT /|| L 2, where a > 0. So that our theory applies to a class of general systems (involving, 
for example, vector equations with matrix operators), we assume the existence of self-adjoint 
"weights", W- and w + such that 

(i) w + is defined on a dense subspace of Tt and on which w + > cl, c > 0. 

(ii) u>_ G OH) 

(iii) u>_ w + P c = P c = P c u>_ w + 

In the scalar case w + and w_ correspond to multiplication by (x) u and (x)~ a , respectively; 
see section 3. 
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The following hypothesis and Corollary ensure that the equation satisfies sufficiently strong 
dispersive time-decay estimates. 

(H3) Local decay estimates on e~ lHot : Let r > 2 + s, where e > 0. There exist w + and W-, as 
above, such that for all /: 

(a) \\w_e- iHot PJ\\ < C (t)- r+1 \\w + fl (2.3) 

and 

(b) IKe-^^o-Ao-^-iO)- 1 P c /|| < C(t)- r+1 \\w+f\\. (2.4) 

(1) Typically one is not "handed" the singular estimate ( |2.4| ). It can however be proved using 
( |2.3|) and the following: 

(bA) Let A denote an open interval containing Ao + fi and contained in the interior of the 
continuous spectrum of H . Let gA denote a smooth function of compact support which is 
identically one on A. Then, for all / such that w + f £ H. 

Ik-e-^'gA /II <C A (t)- r \\w + f\\. (2.5) 



(2) To prove that (|2.4|) is implied by estimates ( ]2.3|) and (|2.5|) , one can follow very closely the 
approach taken in ||37|| , ppj. In particular, see Proposition 2.1 and Appendix A of ||38|| . In [58], 
the proofs are carried out in for a particular choice of TC, w + , and W-, but as indicated in the 
remark following the main theorem of [5_8j , the proof can be adapted to a more general setting 
and is virtually unchanged; one need only replace L 2 by a general Hilbert space 7i, (x)~ a by 
w_, and (x) a by w + . This requires hypothesis: 
(H4) There is a choice of real number c such that 

IKtFo + c)- 1 ^!!,^) (2.6) 

can be made sufficiently small. 

In the scalar case w + and W- correspond to multiplication by (x) cr and (x)~ a , respectively; see 
section 3. 

(H5) Hypotheses on the perturbation W(t): 

We consider time-dependent symmetric perturbations of the form 

W{t) = cos /it /3, with /3* = (3. (2.7) 

In many applications, (5, is a spatially localized scalar or matrix function. 
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To measure the size of the perturbation W, we introduce the norm 

lll^lll = {\\w + P\\ a £M + {\\w+ Pw+Wlw)* ■ ( 2 - 8 ) 

Our final hypothesis is a resonance condition which says that /i + Ao G Ccont(-ffo) and that 
there is nontrivial coupling. 

(H6) Resonance condition There exists 9 > 0, independent of W, such that 

r = - 5(H -X -fi)^ )>6 \\\W\\\ 2 >0 (2.9) 



Remark: That the expression in ( |2.9| ) and the principal value term in ( |2.13| ) below are finite is 



a consequence of the local decay estimate (H3) and the representation: 

( /fyo, {Ho-X-fi- tO)- 1 Pi[> ) = z jH ( Pi/*, e-* 110 -^ 11 ^ ) ds. (2.10) 

Our main result is the following: 

Theorem 2.1. Let H satisfy hypotheses (H1)-(H6) and assume w + <po G TC. Then, if \ \\W\\\ 
is sufficiently small, the solution of satisfies the local decay estimate: 

\\w- 4>(t)\\ < C(t)- r+1 \\w + (j) \\, t G JR. (2.11) 

Remarks: In the proof we actually use a less stringent property than (H6). Rather, we use 
that: 

(H6') III^IH 3 T" 1 is sufficiently small. 

Under the same hypotheses as Theorem 2.1, we obtain more detailed information on the 
behavior of <f>(t): 



Theorem 2.2. 



(f>{x,t) = a(t)ip + <pd{t), (V'o , <Pd{t)) = 0, 
a{t) = a(0) e~ rt e~ tuJl{t) + R a {t) 
Mt) = e~ iHot P C O + 4>(t). 



For any fixed Tq > 0, we have 



\R a (t)\ < C \\\W\\\, |t| < Tor- 1 (2.12) 
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Moreover, 



\R a (t)\ = 0((ty r+1 ), |t|-oo 
LUtit) = \ t - -(/3Vo,(^o-A + ^) -1 Pc^o) t 



4 

- i (/fyb, P.V.(tf - A - v)- 1 PcWo) * + M(t), (2.13) 

wiiere M(t) is of order |||W||| 2 and is oscillatory in t. Furthermore, <f) = <pi + <p2 J 's explicitly 
given in ( |4.S| ), with \\w-(j)(t)\\ = 0((t)~ r+1 ) as \t\ — > oo, and therefore, by (H3) same holds 
\\w-<f> d (t)\\ = 0((t)- r+1 ) as \t\ -> oo. 

Remark: Suppose the initial data is given by the bound state of the unperturbed problem, 
i.e. <fi(x,0) = ipo(x), a(0) = 1, 0^(0) = 0. Then, from the expansion of the solution we have 
that for < t < T^ 1 that the projection of the solution on ip is of order e~ rt , with an error of 
order Hence it is natural to view the state ip e~' lXot as a metastable state, with lifetime 

r = r- 1 ~ \\\W\ 
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3. An example: the Schrodinger equation 

In this section we apply the results of Theorems 2.1 and 2.2 to the special case of the Schrodinger 
equation with a time-periodic and spatially localized perturbing potential, equation ( |1.3| ). Mod- 



els of the sort considered in this example occur in the study of ionization of an atom by a 
time- varying electric field; see fl9[| , 0. 



We begin by indicating assumptions on V(x) and (3{x) under which (H1)-(H6) hold. 

We take TC = L 2 (IR"), and H = —A + V(x), where V(x) and (3{x) are real-valued, 
sufficiently differentiable and decaying rapidly as x G IR n tends to infinity. Thus, Hq is self- 
adjoint, and Hq and W(t) are densely defined operators in L 2 . Moreover, we assume that Ho 
has exactly one negative eigenvalue and the absolutely continuous spectrum of H is equal to 
the positive half-line. 

With regard hypothesis (H3), we take the weights used to measure local energy decay to 
be w± = (x) ±cr . There is a good deal of literature on local energy decay estimates for e~* Hot P c . 
These results require sufficient regularity and decay of the potential V(x). We refer the reader 



to O, pi and p3f; see also [0], M\ 



The large time expansion of the operator: (x)~° e~ ot P c (x)^ cr in L 2 contains a non-generic 
leading order term of order t~z. More rapid decay as t — * oo is ensured provided we assume the 
following condition concerning the behavior of the resolvent of H Q near zero energy: 
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(Z) Zero is not a resonance of the operator —A + V(x). For a precise definition of zero 
energy resonance in this context and a detailed discussion see |U| and |23| . 



Remark: The condition (Z) holds generically in the space of potentials. 
Hypothesis (H3a) is a consequence of: 

Proposition 3.1. Assume that V(x) is sufficiently regular and decays sufficiently rapidly as 
\x\ — > oo; see |T5|/ , ^3jj . Additionally, assume condition (Z). Then, we have the following 
local energy decay estimate. 

(a) Assume n > 3 and o > 1. Then, there exists r > 2 such that 

|| (a;)"' e-^ ot P c / || 2 < C (t)- r+1 \\( X yf\\ 2 . (3.1) 



(b) Assume n = 1. Then, we get ( \3.1\ ) with (t) r+1 — (t) 2. 

(c) Assume n = 2. Then, we get ( |3. 1\ ) with (t)~ r+1 replaced by (£log 2 i) -1 . 



Hypothesis (H3bA) is a consequence of: 

Proposition 3.2. Assume tie above hypotheses on H and assume that A denotes an open and 
bounded subinterval of the continuous spectrum of H . Let #a denote a smooth characteristic 
function, as defined in section 2. Then, for 2 < r < a , 

IK*}- e-^g A {H )f || 2 < C (t)- r \\(x)°f\\ 2 . (3.2) 

proof: If Hq is replaced by the Laplacian, this estimate is straightforward since g/\{Ho)f is 
spectrally supported away from points of stationary phase. An approach to this estimate for 
H , a variable coefficient operator is reviewed in Appendix A and is based on the Mourre 
estimate, where the operator A can be taken to be: A = (x ■ p + p ■ x)/2, where p = —iV. 
We take the time-dependent perturbation to be of the form: 

W(x,t) = cos fit p(x), (3.3) 

with (3 sufficiently differentiable and rapidly decaying in x, e.g. {x) 2rT /3(x) G L 2 . Thus, (H5) is 
satisfied as well. Finally, the resonance condition (H6), holds generically in the potential V(x). 

Therefore, our main theorems Theorems 2.1 and Theorems 2.2 on the structural instability 
of the unperturbed bound state, and large time behavior of such systems apply. 
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4. Decomposition and derivation of the dispersive normal form 



As in [[IT], [P5|l , we begin by deriving a decomposition of the solution, <f>(t) , which will facilitate 
the study of its large time behavior. Let 



Note therefore that (fid = P c 0d since H has been assumed to have only one eigenvalue. 
Remark: The required modifications in our analysis for the case where W(t) has more general 
periodic or quasi-periodic time dependence, and when H has more than one bound state can be 



a sum over the M discrete modes with amplitudes cij(t), j = 0, M — 1, plus the continuous 
spectral part. If W{t) has more general periodic or quasi-periodic time dependence, in the 
analysis of the coupled equations for and (fid{t) ( see below), we use an expansion of W{t) 
in terms of operators of the form (|2.7|) . 

We proceed by first inserting (f4.1| ) into (|2.1|) , which yields the equation: 



(Pit) = a(t)ip + <j) d (t), 



(4.1) 



with the orthogonality condition 



(ipoAd) = for all t. 



(4.2) 



seen at this point. If H has M discrete eigenvalues, then the decomposition ( f4.1| ) is replaced by 



id t a(t)ip + id t (f>d{t) 



= X a(t)^ + Ho4> d (t) 

+ a(t)W(t)if> + W(t)Mt) 



(4.3) 



Taking the inner product of ( |4.3| ) with ip yields the following equation for a(t): 




X a(t) + (i/) ,W(t)i/) )a(t) + i$v,W(t)<t> d ), 
(^o,0(O)) 



(4.4) 



In deriving ( |4.4|) we have used that ip$ is normalized and the relation 



(V> , d t (j) d ) = 0, 



(4.5) 



a consequence of ( |4.2| ). 

Applying P c to ( |4.3| ), we obtain an equation for (fid'- 



id t (fid{t) 
MO) 



H Mt) + P c W(t)Mt) + a(t)P c W(t)ifi , 
P C 0(O) 



(4.6) 
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Since we are after a slow resonant decay phenomenon, it will prove advantageous to extract the 
fast oscillatory behavior of a(t). We therefore define: 

A(t) = e tXot a(t). (4.7) 



Then, ( |4.4| ) reads 

d t A = -iA(iJ> ,W(t)ipo) - te iXot (^W(t)Mt))- (4- 
Solving (^4.6|) by Duhamel's formula we have 

Mt) = e- iHot <f) d (0) - i C e- iHo{t - s) V c W{s)a{s)^ds 



- i J e- iHo(t - s) P c W{s)<j) d {s) ds 

= Mt) + Mt) + Mt)- (4.9) 

By standard methods, the system (gg)-(gg) for A(t) and Mt) = M) ~ e ~* Ao< Mt) has a 
global solution in t with 

A e C\M), \\Mt)\\ e C°(1R), \\w-<j> d {b)\\ e C°(1R). 

Our analysis of the \t\ — > oo behaviour is based on a study of this system. 
We next insert (|4.9| ) into ft4.8|) we get 

d t A(t) = -tA(t)^ ,W(t)M - te lXot j2(MW(t)M- (4-10) 

j=0 

We next give a detailed expansion of the sum in Qll(J| ). It is in the j = 1 term that the key 
resonance is found. This makes it possible to find a normal form for ( |4.1(J| ) in which internal 
damping in the system is made explicit. This damping is responsible for the energy transfer from 
the discrete to continuum modes of the system and the associated radiative decay of solutions. 

Proposition 4.1. 

-ie iXot cos fit ( fiipo, Mt) ) = R1 + R2 + R3 + R±, where 
R 1 = ~ e iXot cos fit J* (j3ip , e - iHo{t ' s) e i{p - Xo)s A(s)V c (3%) ds 



R 2 = l -e' iflt cos fit A(t) ^^(Ho-Xo-fi-iO^PM 
• ^ J*ot nno ,,± ! 0.1. c -iH t(TT -v _ -n^-ii 



R 3 = --A(0)e iXot cos fit [Mo, e^ ot (# - A - /x - iOpPc/fyt 

R 4 = - l - e iXot cos fit J' (A/>o, e- iHo ^~ s) {H - A - fi - i0Me- i{Xa+ ^ s d s A{s)P c (3^ ) ds. 

(4.H) 
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Proof of Proposition 4.1 

Using that cos at = |(e* M * + e~ tM *) and the definition A(t) = e lXot a(t), we get from ( |4.9| ) 

<f>i(t)= -I f* e- iHo{t ~ s) e- i{x °-' M)s A{s)P c ^ ds 
Jo 

-f f e- iHo{t - s) e- i{Xo+,l)s A(s)P c pip ds 

= K + K ves . (4.12) 

Therefore, 

-ie iXot cos at (J3if)o,<f>i(t)) = -ie iXot cos at (J3ifo, K + K ICS ) 

= R x - ie iXot cosa.t(pip ,K ms ) (4.13) 



The terms Rj, j = 2,3,4 in ( |4.11 ) come from an expansion of the last term in (|4.13|) . This is 



the source of the key resonant contribution. We now proceed with a careful evaluation of this 
term. 

We first regularize K res . For t > 0, let e be positive and arbitrary and define: 

K £ res = - l - f e - iHo ^e- l(Xo+ ^- i£)s A(s) P c Mo ds (4.14) 



2 



o 



Note that K ies = lim £ ^ -^res- Integration of the expression in ( [4. 14| ) by parts, and letting e tend 
to zero from above gives 

Proposition 4.2. The following expansion of K res is valid in S': 
K Tes = -I e -^o+M)^( t )( Fo _ A _ / ,_ i0 )-i Pc ^ 
+ ±A(0) e- iHot (H — X — fi — 20)- 1 P c /?Vo 

+ I f e- iHo(t - s \H - \- u.-iQ)- 1 e- l{Xo+ < l)s d s A(s)P c Ptp ds (4.15) 
2 jo 



The choice of regularization, — ie, in (|4.14 ) ensures that the latter two terms in the expansion 



of K Ies , (|4.15|) , decay dispersively as t — > +oo; see hypothesis (H3) and section 6. For t < 0, 
we replace +ie with —is in (|4.14|) . 



The proof of Proposition 4.1 is now completed by substitution of the expansion Q4.15 ) for 
K res into the last term in ( |4.13|) . 

We further expand the terms Rx and R2 to make explicit the resonant contributions. 
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Proposition 4.3. 



R x = p%(t)A{t) + R la (t) + R u (t), 
R 2 = (p 2 (t) - T)A(t), 



(4.16) 
(4.17) 



where 



P2(t) 

r 

Rla{t) 
Rlb{t) 



-(1 + e 2 ^ 
4 



:i + e" 



2 ""W w- ,P.V. (//„ A p) 1 P c ^o) 



7T 

- 4 6 

7T 



-2i/it 



(Mo,5(H -A Q -//)P c /3Vo) 



(/fyo,*(#o-A -//)Pc/tyo) 



-^(0)e iA °'cos/xt e-^ ' (F - A + /i^P^o) 



2i 



-e tXot cos, 



« /it (Wo, jf e-^ o{ *- s) e- i(A °-^ s 9 s A(s)(iJo - A + p)~ l P c (3^ 
Proof. To obtain ( [4.17] ) we use the well known distributional identity 

(x - iQ)- 1 = py. x- 1 + m5(x). 



Equation (|4.16|) follows by integration by parts. 

Combining equation ( |4.10| ) with the previous propositions, we have 



Proposition 4.4. 



d t A(t) 



-r + pit) ) A{t) + E(t), 



where 



E(t) = R la + R u + R 3 + R 4 - ie iXot cos pt M*)) ~ ^ iA °' cos pt (/^ , 2 ) . (4.20) 



iXot 



Here, O an d 4>2 are given in (\4.9j ), Ri a and Ru in ( \4.18j) , R3 and R4 in ( \4.1J\ ), and 



pit) = -i cos pt (ipo, /3ipo) + p 1 (t)+p 2 (t) 
= ^(/3^o,(^o-A + /i)- 1 Pc/3^o 
+ ^(/3Vo,P.V. (# - A - /i^Pc/^o) + m(t), 



(4.18) 



(4.19) 



(4.21) 



where m{t) is a bounded and oscillatory function of zero mean. 



Remarks: 

(1) The point of ( 4.19|) is that the damping coefficient (T > 0), which arises due to the coupling of 



the discrete bound state to the continuum modes by the periodic perturbation, is made explicit. 
This facilitates a direct study of the transfer of energy from the bound state to the continuum 
and the associated radiative decay of solutions. 

(2) The leading order part of equation ( |4.19| ) is the analogue of the dispersive normal form 
derived in |33 for a class of nonlinear dispersive wave equations. 



In the next section we estimate the remainder terms in (|4.9| ) and ( }4.19|) . 
5. ODE Estimates 



Equations ( [4.9| ) and ( |4.19[ ) comprise a dynamical system governing 4>d(t) and a(t), the solution 



of which is equivalent to the original equation (1.1). Our aim in this decomposition is to 
separate the dispersive (polynomially time-decaying) part 4>d{t) from the solution, as well as the 
transient (exponentially time-decaying) part. In this and in the following section we derive a 
coupled system of estimates for A(t) and 4>d(t), we show that A{t) decays in time, provided 4>d(t) 
is dispersively decaying and vice-versa. Then, we exploit the assumption that the potential, W 
is small to close the resulting inequalities, and prove the decay of both A(t) and 0d(t). 

At this point, we state a simple Lemma which we shall use in a number of places in this and 
in the next section. 



Lemma 5.1. Let a > 1. 

t-t 



o 



(t - S)~° (S)-? ds < C a> p (t)-min(ft/3+ Q -l) (5.1) 



Proof: The bound is obtained by viewing the integral as decomposed into a part over [0, t/2] and 
the part over [t/2, t]. We estimate the integral over [0, t/2] by bounding (t — s)~ a by its value 
at t/2 and explicitly computing the remaining integral. The integral over [t/2, t] is computed 
by bounding (s)~@ by its value at t/2 and again computing explicitly the remaining integral. 
Putting the two estimates together yields the lemma. 

We now turn to the estimate for A(t) in terms of the dispersive norm of <pd{t) and local decay 
estimates for e~ lHot P C (H ) . 



Proposition 5.1. A(t), the solution of ( |4.1fy ), can he expanded as: 



A(t) = eti pis)ds (e- rt A{0)+R A (tj) (5.2) 
R A (t) = f e~ r{t - T) E(t) dr, (5.3) 
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where E(t) is given in ( <L2G ) and ( |5.Sp . For any a > 1, R A (t) satisfies the estimates: 



sup (t) 1 - 1 \R A (t)\ < C ie - l ~ sup \E(t)\ + CT' 1 sup (r^l^r)! ) ,(5.4) 
2F' a <t<T o<r<r- a r- a <r<T 

sup (t) 7 - 1 \R A {t)\ < DT- ar sup |£(r)| (5.5) 

o<t<2r- a o<r<2r- a 



(5.6) 



Proof. To prove ( |5.4p we begin with (4.19) 

d t A(t) = p(t)A(t)-VA(t)+E(t), 
where E(t) is given by ( |4.20[ ). Let 

A(t)=e~Io p ^ ds A(t). 



Then, A satisfies the equation 



d t A = —TA + E(t) 
E(t) = e-ti p{s)ds E(t). 



Solving; (15.81) we get 



A(t) 



e- rt A(0) + C e- Tit - s ^E(s) ds 
Jo 



= e- rt A(0) + R A (t). 



(5.7) 



(5.8) 
(5.9) 



(5.10) 
(5.11) 



Since the real part of pit) is of the general form of C\ coscuit+c^ sinco> 2 £, its integral is bounded 
uniformly in t. Therefore, 

C^li^)! < \A{t)\ < C\A(t)\ for some C > 0, 

and similarly for E(t),E(t). It is therefore sufficient to estimate A(t), in terms of E(t). 
Remark: Estimates of R a {t), which appears in the statement of Theorem 2.2, are related to 
those for R A {t) via : 



\Ra(t)\ 



f p(s) ds R ^ 



< C \R A {t)\ 



(5.12) 



Note that on any fixed time interval, < t < 2M, one has by Gronwall type estimates that, 
\A(t)\ and ||w_0d(t)|| are bounded by a constant depending on the initial data and M. This 
bound is, in general, exponentially large in M. We shall now focus on obtaining appropriate 
bounds for these quantities on t— intervals, [2M, T), which are independent of T. 
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From (|5.10|) we have for any M > 0: 

\A(t)\ < \A(0)\e- rt + e- r ^ s) \E(s)\ds+ e - r(t - s) \E(s)\ ds 

JO JM 

= \A(0)\e- Tt + h(t) + I 2 {t). (5.13) 

Set 

M = T~ a , a>l. 
We now estimate the terms hit) and hit) in ( |5.13| ) for 2T~ a <t <T. 

( t y- 1 hit) = (ty- 1 / e^^lEis^ds 

Jo 

< (t) r - 1 e-3 r * . f M e' r ^- s) ds ■ sup \E(r)\ 

Jo o<T<r- a 

< sup i(t) r - l e- r i) ■ CT' 1 • sup |%)| 

2T- a <r<T 0<T<T- a 

< Ce~ r ~ S sup \E(r)\, (5.14) 

o<r<r- Q 

for some 5 > 0. Therefore, 



sup {{ty- 1 hit)) < Ce- L sup \E(t)\ (5.15) 

2T- a <T<T 0<T<r~ a 

We estimate hit) on the interval 2T~ a < t < T as follows: 

(tY- 1 h{t) < (ty- 1 f e-^ l - s \s)- r+1 ds sup Ur) r - l E{r)) (5.16) 

Jr- a r~ a <T<T 

The integral is now bounded above using the estimate 

(ty- 1 f e~ r{t - s) (s)- r+1 ds < CT-\ t > 2T~ a . (5.17) 

This gives 

sup (t)^ 1 hit) < CT' 1 sup ((ry^Eir)) (5.18) 
2T- a <t<T r~ a <r<T 



Assembling the estimates (|5.15| ) and ( 5.18|) yields estimate ( |5.4|) of Proposition 5.1. Estimate 



3|) is a simple consequence of the definition of Ra{P)- 
6. Dispersive Estimates and Local Decay. 

In this section we prove the local decay of 0<j and the decay in time of the remainder terms in 
the ODE of Section 4. To this end we will repeatedly use the following: 
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Lemma 6.1. For any rj G [0, r — 1] 

[ l w_e- iHo{t - s) PJ(s)ds 
Jo 



< C(t)-« sup «t>"|K/(t) 

0<T<t 



and 



f w.e- m °^P c (H - A - n - i0)- x /(s) ^ 

JO 



< C(t>-" sup «t>" |K/(r) 

0<T<t 



(6-1) 



(6.2) 



Proof. The proof follows from the assumed local decay estimates on e see (H3). Namely, 
using that r > 2, 

J' w^e- iHo{t - s) PJ(s) ds < j l \\w.e- iHo{t - s) P c W-\\ c{H )(s}' v ds 

■ sup «t)"|K/(t)||) 

0<t<* 

< C [\t-s}- r+1 (s)-*> ds sup ((r)"|K/(r)||) 

< sup ((rr||«; + /(r)||) 

0<r<< 



which proves (fTT|). The proof of ( |6.2| ) is identical, and uses the singular local decay estimate of 
(H3). 

We now define the norms 

(6.3) 



and 



[A] a (T) = sup (rr\A(r) 

0<t<T 



[4>d\LD,a{ T ) = SU P ( r ) a H W -0d( r ) 
0<t<T 



(6.4) 



Then we have 

Proposition 6.1. For any T > 0, 

[Mld,t-i(T) <C( \\w + <f> d (0)\\ + III^IH [A] r _i(T) ). (6.5) 
Proof. From equation (|4.9|) we get, using the assumed local decay estimate for e~ lHot and (|6.1|) , 



\v>-Mt)\\ < E 

j=0 

< C(t)- r+1 \\w + U0)\\+ C{t)- r+1 [A] r _i(*) sup 



0<s<t 



+ c infill (t)- r+1 [M L n M (t). 



(6.6) 
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Since ||w + W / (s)-?/'o|| < 1 1 \W\ \ \ \\ipo\\ = \\\W\\\ and || \W\ || is assumed to be small, multiplying 
both sides of this last equation by (t) 1 "" 1 and taking supremum over t <T yields ( |6.5|) . j 
We now estimate E(t). 

Proposition 6.2. Let T > 0. For any r] G [0, r — 1]: 

[£WT) < c( infill 2 |A(o)| + iiiiyiiin^wii + infill 3 [A],(r)) (6.7) 

Proof. 



E(t) is defined in ( |5.6|) and ( |4.19j ). From these equations it is seen that we need to estimate 
the following terms: 



(aij 

(a 2 ) 

(«s) 

(a 4 ) 
(«s) 



fl la | = \A(0)\ (f3ip ,e- iHot (H -X + fi)- 1 P c ^ 



R 
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1^4(0)1 
1 



o,e 



-iH (t-s) i(n-\ )s 



(H -\ + t i)- 1 d s A(s)P c P^ ) ds 



R 3 \ = -\A(0)\ (/3Vq, e-^ *(F - A - /i - zO)^P c ^ ) 



Wo, e 



-iH (t-s) 



Wo,0o)| = |Wo,e-^Vd(0))| 



(F - A - // - 20)- 1 e- i ( Ao+ ^ s 9 s A( S )P c /3Vo) 



(^0, 02^)1 



(/?Vo, e" 



-iH (t-s) 



cos /xs P c P4>d(s) ds) 



The estimates of (a.,) repeatedly use Lemma 6.1. Let r] G [0, r — 1]. 
Estimation of («i): 



< |-4(0)| (w+^Vo, W-e~ iHot (H - A + ^-'Pcw- w+Pfa) 

< c\a(o)\ \\\w\\\ 2 {ty r+1 



(6.8) 



by local decay estimates and our assumptions on H , since Ao — p is not in the spectrum of H . 
Estimation of (a 2 ) 



From (15.61) we have that 



\d s A(s)\<C\\\W\\\ \A(s)\ + \E(s) 



(6.9) 



since p and V are, respectively, linear and quadratic in 
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Applying Lemma 6.1 to (0J2) we then get 



\Rib\ < |-4(o)| 



(w + (3ip , w_e- lHo{t - s) (H - A + / u)- 1 P c w_ d s A(s) w + (3ip ) ds 



< C\\\W\\\ 2 (t)-m\\W\\\ [A] v (t) + [E] n (t)). (6.10) 

Estimation of (a 3 ) Since Ao + /x G ct(Hq), we use (|2.4|) together with the procedure used to 
estimate (a 2 )- The result is: 



\R 3 \ <C\A(0)\ \\\W\\\ 2 (t)- r+l . (6.11) 

Estimation of (0*4) 

This term is estimated in a manner similar to 02 and oi%: 

\Ra\ < c\\\W\\\ 2 {\\\W\\\ [A} v (t) + [E] v (t)} (*)-". (6.12) 

Estimation of (0:5): 

Since, by definition, <f) a (0) = P c 0d(O) we can apply local decay estimates for e~ tHot to get 

m ,Mt))\ < C\\\W\\\ (t)~ r+1 \\w + M0)l (6.13) 

Estimation of (a^) 

Applying Lemma 6.1 as before we get, for < t < T, 

|(/tyb,0a(f))| < C\\\W\\\ 2 (t)-*[<j> d ]LD, v (T). (6.14) 

Using Proposition 6.1 to estimate \4>d\LD,r,{i) in ( |6.14|) , we get 

m^Ht))\<C\\\W\\\ 2 (t)-*{ \\w + U0)\\ + \\\W\\\ [A] v (t) }. (6.15) 

Finally, combining the above estimates for (acj), 1 < j < 6, we can now bound [E^T) for 
any 77 G [0, r — 1] as follows: 

[E] v (t) < C [\\\W\\\ 2 |A(0)| + III^IH ||w+^(0)|| + III^IH 2 [E] V (T) + III^IH 3 [A] n (T)} . 

(6.16) 

Since |||W||| assumed to be small, Proposition 6.2 follows. [| 

We can now complete the proofs of our main results, Theorem 2.1 and Theorem 2.2. To 
prove the assertions concerning the infinite time behavior, the key is to establish local decay 
of 4>d, in particular, the uniform boundedness of [4>d]LD,r-i(T). This will follow directly from 
Proposition 6.1 if we prove the uniform boundedness [A] r _i(T), or equivalently [A] r _i(T). 
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Proposition 6.3. There exists a constant C*, depending on ||0o||, ||w+0o||, 1 1 \W\ | |, V and r, such 
that for any T > 

[A] r -i(T) < (6.17) 

Proof. 

We begin with the expansion of A(t) given in Proposition 5.1. Multiplying ( |5.2j ) by (t) r ~ 1 , 
and taking the supremum over < t < T we have: 

[A] r -i(T) < (\A(0)\ T~ r+1 + sup \R a (t)\ + sup \R A (r)\ ) (6.18) 

\ 0<r<2r- a 2T' a <T<T I 



The right hand side of (|6.18|) is estimated using Proposition 5.1. 



[A] T - X (T) < (|A(0)| + D T- ar [E] (2T- a ) + C x e~ r ~ S [E] (2T- a ) + C T-^E^T) 

(6.19) 

Next, we apply Proposition 6.2 which yields: 

[A] r -i(T) < (\A(0)\ T- r+1 + D T- ar [E) (2T- a ) + [E] (2T- a )) 

+ CT- 1 (\A(0)\\\\W\\\ 2 + \\\W\\\ \\w + MV)\\ + \\\W\\\ 3 [A] r -i(T)) . 

Note that by Proposition 6.2 and the simple bound: 

[A] r -i(T) < C (T)'- 1 [A) (T) < C (Ty- 1 ||0 O ||, (6.20) 

[i?]o(2r _a ) is bounded in terms of the initial data and |||W|||. 
By hypothesis (H6), there is a positive constant 9q, such that 

iiim^iii 3 iiiwiii 

INJL < NUN. (6 . 21) 

Therefore, for |||W||| sufficiently small 

[A] r -i(T) < C, (6.22) 



(Alternatively, we can obtain ( 6.22j) by imposing that |||W||| 3 r 1 be sufficiently small.) Here, 



depends on ||0o||, ll^+^oll, lll^lll; r and r. 

This completes the proof of Proposition 6.3 and therewith the t — > oo asymptotics asserted 
in Theorems 2.1 and 2.2 

Remark: This estimation procedure differs from one used in |3*8" |, where the cases r — 1 large 
and r — 1 small are treated differently. 
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R A {t) <C /"* e- r (*- r ) |£(r)| dr. 
Jo 



(6.23) 



It remains to prove the intermediate time estimate (|2.12 ). The ingredients are contained in 
and its proof. First, by (|5.3| ) 

ft 
10 

Let T denote an arbitrary fixed positive number. We estimate (|6.23|) for t G [0,Tor -1 ]. We 
bound the exponential in the integrand by one (explicit integration would give something of 
order T -1 ), and bound |-E(t)| by estimating the expressions (<x,), j = 1,...,6 in the proof 
of Proposition 6.2. First, the estimates of Proposition 6.2 for ai,a^ and a§ are useful as is. 
Integration of the bounds ( |6.8| ), ( |6.11D and ( |6.13|) gives: 

I, j = 1,3,5. (6.24) 



-r(t-r) 



>i) < C\\\W\ 



To estimate the contributions of (a 2 ) and (0:4), first observe that by (|6.9|) and Proposition 6.2 
with f] = 

\d s A(s)\ < C \\\W\\\ ||w + 0(O)|| 



(6.25) 



Therefore, using local decay estimates we have: 

/ e~ r{t ~ T) (a,-) dr < C ToT' 1 \\\W\\\ 3 \\w + <f>(0)\\ 
Jo 



< 



\W\ 



3 = 2,4. 



Finally, we come to the contribution of 



"6j 



We rewrite (a 6 ) as follows. 



((3 e lHo(i - s) P c /3Vo,cos^ <j) d (s) ds) 

j {w + j3w + ■ w^e iH °^~ s ^P c w^ ■ w + f3ijj 0} cos /js w_0^(s) ds) 
Jo 



(6.26) 



Recall that by ( f4.9| ) <f) d = <f) + (pi + <p 2 , where <po(t) = e~ zHot <pd(0). Using local decay estimates 
(H3), the contribution of the term (po(t) to (a 6 ) can be bounded by C \ \ \W\ | | 2 Hw+^O) || (r) _r+1 . 
Multiplication of this bound by e~ r ^~ T ^ and integration with respect to t gives the bound 
C \ \\W\\\ 2 \\w + (()d(Q)\\- To assess the contributions from <pi + 2 , note that local decay estimates 
(H3) imply 



||W-(01+0 2 )|| < C \\\W\\\ \\w + (p \\. 
Putting together the contributions from (po and from (pi + (p 2 , we have: 



-r(t-r) 



[a 6 ) dr < C 



\W\ 



\w + <PM\\ 



v- 1 \\\w\ 



(6.27) 



(6.28) 



The above estimates and (H6) imply ( [2 . 1 2| ) . This concludes the proof of Theorem 2.2. 
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7. Appendix A: General approach to local decay estimates 



Hypothesis (H3) for our main theorem is one requiring that our unperturbed operator, Hq, 
satisfy a suitable local decay estimates. It is remarked following this hypothesis that, in practice, 
the verification of (|2.4j) using (|2.5| ). In this section we give an outline to a very general approach 
to obtaining estimates of this type, based on a technique originating in the work of Mourre 



see also ||25|| . In the following general discussion we shall let H denote self-adjoint operator on 
a Hilbert space, 7i, keeping in mind that our application is to the unperturbed operator H . 
Let E e o~(H), and assume that an operator A can be found such that A is self-adjoint and 
T>(A) PI Ti is dense in 7i. Let A denote an open interval with compact closure. We shall use the 
notation: 

&d n A (H) = [...[H,A],A],--A], (7.1) 

for the n-fold commutator. 

Assume the two conditions: 
(Ml) The operators 

g A {H) [H, A] , [H, A] g A (H), and g A (H) ad™(#) g A (H), 2<n<N (7.2) 
can all be extended to a bounded operator on 7i. 
(M2) Mourre estimate: 

g A (H) i[H, A] g A (H) > 9 g A (H) 2 + K (7.3) 
for some 9 > and compact operator, K. 

Remark: The theorems below were orginally proved under more restrictive conditions than 



(Ml) H, M, p5|, p4|, |36|. The more general results stated here can be proved using the 



approach to velocity bounds in [12 . 



Theorem 7.1. (Mourre; see ||/ Theorem 4.9, Lemma 5.4) 

Assume conditions (M1)-(M2), with N=2. Then, in the interval A, H can only have absolutely 
continuous spectrum with finitely many eigenvalues of finite multiplicity. Moreover, the operator 

(A)- l 9A (H) (H-z)- 1 (A)- 1 (7.4) 

is uniformly bounded in z, as an operator on 7i. If K = 0, then there are no eigenvalues in the 
interval A. 
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Theorem 7.2. (Sigal-Soffer; see |3§ £2}) 

Assume conditions (M1)-(M2) with N > 2 and K = 0. Then, for all e > 

|| f(j<9) e~ lHt g A (H)^ || 2 < C(t)- N+1+e \\ (A)""V || 2 , (7.5) 

and therefore 

|| (A)~ a e~ iHt g A (HW || 2 < C (t)- \\ (A) N ~ 1 ip || 2 , (7.6) 

for a < N — 1. Here, F is a smoothed out characteristic function, and F (j < is defined by 
the spectral theorem. 



Remark: a is required to be smaller than y in H, pi and smaller than that in [B(|. The 



above bound a < N — 1 is proved in [ 121 ; see Theorem 1.1 and the remark below it. 



Let Ai denote an open interval containing the closure of A. 

Corollary 7.1. Assume that (x)~ a g Al (H)(A) a is bounded. Then, in the above theorems we 
can replace the weight (A)~ a by (x)^ cr . 

The strategy for using the above results to prove local decay estimates like that in (H4) is 
as follows. Then 

|| (x)^ e~ lHt g A {H)^ || 2 = || (x)-"g Al (H)e- iHt g A (H)il>\\ 2 

= || (x)-°g Al (H)(Ay ■ (A)- e~ im g A (H)iP || 

< || (x}^g Al (H)(Ay\\ ■ || (A}~°~ e- im g A {H)^ \\ 

< C || (A)- e- iHt g A {H)^ \\ 

< d|| F < 0) (A'j u e~ lHt g A (H)^ || 2 



t 

+ C 2 || F (y > 9 ) (A)-° e- im g A {H)^ || 2 (7.7) 

Theorem 7.2 is used to obtain the decay of the first term on the right hand side of ( |7.7| ), while 
we can replace A by 6t in the second term. 
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